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Foreword 

The  calculations  in  this  note  have  a  form  similar  to  those 
in  a  few  previous  notes  on  cylindrical  loops.  This  note  extends 
these  types  of  calculations  to  spherical  geometries  where  the 
sensor  in  this  case  measures  the  displacement  current  density. 

For  convenience  the  figures  are  grouped  together  after  the  sum¬ 
mary  and  before  the  appendices.  Appendix  C  was  written  by  Mr. 

Joe  P.  Martinez  of  Dikewood  and  we  would  like  to  thank  him  and 
Mr.  Larry  Berg  of  AFWL  for  the  numerical  calculations  and  graphs. 


Contents 


Section  Page 


I.  Introduction  2 

II.  Electromagnetic  Fields  in  Spherical  Coordinates  4 

III.  Vector  Plane  Waves  in  Spherical  Coordinates  11 

IV.  Short  Circuit  Current  19 

V.  Admittances  26 

VI.  Frequency  Response  Characteristics  35 

VII.  Summary  37 

Figures  38 

Appendices 

A.  Calculation  of  An  52 

B.  Properties  of  Admittance  Series  for  Large  n  58 

C.  Numerical  Techniques  for  Computer  Calculations  66 


I .  Introduction 

Among  the  problems  of  electromagnetic  sensor  design  for  use 
in  non-conducting  media  there  are  the  general  problems  of  sensor 
accuracy,  directional  sensitivity  at  high  frequencies,  and  maxi¬ 
mizing  the  upper  frequency  response  for  a  given  sensitivity.  Here 
we  are  concerned  with  sensors  for  measuring  some  electromagnetic 
field  component,  or  its  derivative,  with  a  flat  frequency  response 
Over  the  bandwidth  of  interest  so  that  in  measuring  pulsed  fields 
there  is  no  distortion  of  the  waveform,  within  limits  like  the 
risetime.  As  an  example,  a  multi-gap  cylindrical  loop  has  a  well 
calculable  equivalent  area  for  measuring  a  component  of  B,  and  by 
increasing  the  number  of  gaps  the  upper  frequency  response  can  be 
raised  and  the  directional  sensitivity  at  high  frequencies  is  re¬ 
duced  (for  waves  still  propagating  perpendicular  to  the  loop 
axis) . 1 

In  this  note  we  consider  another  kind  of  sensor  which  we 
term  a  hollow  spherical  dipole.  Since  this  sensor  is  based  on  a 
sphere  we  can  analyze  its  performance  using  vector  eigenfunction 
expansions  in  spherical  coordinates.  (This  in  itself  is  a  good 
reason  for  considering  a  spherical  sensor.)  The  analysis  will 
follow  an  approach  similar  to  that  used  in  two  previous  notes 
concerning  cylindrical  loops  where  cylindrical  vector  eigenfunc¬ 
tions  were  used.1/2 

The  basic  sensor  geometry  is  shown  in  figure  1.  It  consists 
of  a  hollow  sphere  of  radius  a  with  a  gap  of  angular  width  2i(i0 
symmetrically  cut  around  the  equator  of  the  sphere.  The  sensor 
is  described  in  spherical  coordinates  (r,  0,  <j>)  as  a  conducting 
surface  on  r  =  a  for  0  <  0  <  0O  and  ir  -  0O  <  0  £  ir.  The  gap  is 
described  by  r  =  a  and  0O  <  0  <  tt  —  0O  where  0  <  0O  <  ir/2 .  We 
have  the  relation 


(1) 


There  is  also  the  cylindrical  coordinate  system  (¥,  <p ,  z) 
and  the  sensor  geometry  is  constrained  to  have  axial  symmetry 
(about  the  z  axis)  so  that  its  response  is  independent  of  The 
gap  is  assumed  to  have  resistive  loading  uniform  in  <{>  to  preserve 


1.  Capt  Carl  E.  Baum,  Sensor  and  Simulation  Note  41,  The  Multi- 
Gap  Cylindrical  Loop  in  Non-Conducting  Media,  May  1967. 

2.  Capt  Carl  E.  Baum,  Sensor  and  Simulation  Note  30,  The  Single- 
Gap  Cylindrical  Loop  in  Non-Conducting  and  Conducting  Media, 
January  1967. 


axial  symmetry.  This  resistive  loading  might  in  practice  be 
many  cable  inputs  evenly  spaced  around  the  gap;  the  cables  would 
bring  the  total  current  crossing  the  gap  with  equal  delays  to  one 
common  point  where  the  signal  is  desired.  Such  cable  networks 
are  not  considered  in  this  note,  but  they  are  assumed  to  be  lo¬ 
cated  in  positions  which  do  not  significantly  perturb  the  sensor 
geometry . 

The  basic  mode  of  operation  of  this  type  of  sensor  uses  the 
short  circuit  current  across  the  loop  gap.  As  will  appear  in  the 
analysis,  for  wavelengths  much  larger  than  a  the  short  circuit 
current  is  just  3ira2  times  the  z  component  of  the  total  current 
density.  If  the  medium  conductivity  is  zero  then  the  total  cur¬ 
rent  density  is  just  the  displacement  current  density.  (There 
are  no  source  currents  in  these  calculations.)  Thus  one  might 
refer  to  this  sensor  as  a  D  sensor  or  a  total  current  density 
sensor,  depending  on  the  specific  application . 3  The  sensor  has 
an  equivalent  area  of  37ra2  which  is  quite  accurate  as  long  as  'rQ 
is  small.  The  actual  accuracy  of  this  number  for  the  equivalent 
area  is  not  considered  in  this  note.  The  simplifying  assumptions 
allowing  the  present  high-frequency  analysis  give  3iTa2  as  the 
equivalent  area. 

In  outline  this  note  first  considers  the  expansion  of  an  in¬ 
cident  electromagnetic  plane  wave  in  terms  of  spherical  vector 
eigenfunctions.  Then  this  plane  wave  is  imposed  on  the  sphere 
with  a  shorted  gap  in  order  to  calculate  the  short  circuit  cur¬ 
rent  as  a  function  of  frequency  and  the  angle  of  wave  incidence. 
Then  assuming  small  and  a  quasi  static  electric  field  distri¬ 
bution  across  the  gap  we  calculate  the  sensor  admittances  associ¬ 
ated  with  both  the  volume  inside  and  outside  the  sensor.  Combin¬ 
ing  the  admittances  due  to  the  sensor  geometry  and  the  assumed 
cable  loading  with  the  short  circuit  current  then  gives  the  sensor 
response  functions  for  an  incident  plane  wave.  From  these  results 
things  like  optimum  cable  loading  can  be  calculated.  For  this 
sensor  with  the  case  of  interest  that  the  medium  conductivity  is 
zero  the  admittance  due  to  the  cable  loading  is  lar'c  compared  to 
the  sensor  admittance  (basically  a  capacitance)  for  frequencies 
below  the  upper  frequency  response  of  the  sensor.  The  sensor  re¬ 
sponse  is  then  normalized  by  dividing  the  current  output  by  3?Ta2 
times  6,  the  ideal  low-frequency  response. 


3~.  Capt  Carl  E.  Baum,  Sensor  and  Simulation  Note  38,  Parameters 
for  Some  Electrically-Small  Electromagnetic  Sensors,  March  1967. 
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II .  Electromagnetic  Fields  in  Spherical  Coordinates 


Consider  that  we  have  a  linear,  homogeneous,  isotropic  medium 
with  permittivity  e,  permeability  y,  and  conductivity  a.  We  have 
propagation  constants^ 


k  =  /-iwy (o  +  iwe) 

(2) 

y  =  /sy (a  +  se) 
and  a  wave  impedance 

(3> 

where  s  is  the  Laplace  transform  variable  which  we  take  as  io)  for 
the  frequency-domain  analysis  in  this  note.  The  radian  frequency 
is  w  and  i  is  the  unit  imaginary.  Our  interest  lies  for  the  most 
part  in  a  -  0,  as  is  used  for  the  numerical  calculations.  How¬ 
ever  we  keep  a  in  the  analysis  for  generality. 

With  time  harmonic  fields  and  with  ela)t  suppressed  Maxwell's 
equations  have  the  form 

V  *  3  =  -itoS  ,  V  x  h  =  J  +  ico5 

(4) 

V  •  3  =  0  ,  V  *  D  =  p 

together  with  the  constitutive  relations  and  Ohm's  law 

3  =  yH  ,  5  =  ,  3  =  o3  (5) 

From  equations  4  and  5  one  obtains  vector  wave  equations  of  the 
form 

V2E  +  k2E  =  <5  ,  V2H  +  k2H  =  5  (6) 


Note  that  there  are  assumed  to  be  no  source  currents  or  charges 
(p  =  0)  in  the  medium,  but  therg  will  be  charges  (and  currents) 
on  the  spherical  sensor.  Thus  E  (as  well  as  H)  has  zero  divergence 
away  from  boundaries  allowing  the  result  in  equations  6. 


All  units  are  rationalized  MKSA . 


In  spherical  coordinates  the  solution  of  the  scalar  wave 
equation  can  be  written  as  a  linear  combination  of  functions  of 
the  form^ 


m. 


f^£)  (kr)P^(cos  (9 ) ) 


cos  (mcj>) 
sin  (m<j>) 


(7) 


where  f^)  fkr)  is  one  of  the  spherical  Bessel  functions  jn(kr), 
yn(kr),  h^)  (kr)  ,  h^2)  (kr)  for  £  =  l,  2,  3,  4  in  that  order. 

The  third  argument  of  H  is  listed  as  e  or  o  (meaning  even  or  odd) 
e  corresponds  to  using  cos(m<p)  and  o  to  using  sin(m<]>).  Unless 
noted  to  the  contrary  the  definitions  of  the  special  functions 
correspond  to  those  in  a  standard  reference  work. 6  In  particular 
the  Legendre  functions  of  the  first  kind  P™(£)  of  degree  n  and 
order  m  on  the  cut  (-1  <  £  £  1)  in  the  complex  £  plane  have  the 
form  (ref.  6  eqn.  8.6.67 


P™(£) 


<-l)m(l  -  £2) 


dC 


P  (5) 
m  n  ^ 


where 


P  (O 
n 


p°(U  = 


.n 


n 


«n  *  j  n 
2  nl  d£ 


(€  ~  1) 


(8) 


(9) 


where  we  only  consider  n  and  m  as  non  negative  integers  and  £  as 
a  real  argument  with  -1  £  £  £  1.  Our  definition  differs  from 
that  in  some  texts  on  electromagnetic  theory^/7,8  in  which  the 
factor  of  (-l)ra  is  not  included.  The  form  of  equation  8  is  how¬ 
ever  consistent  with  various  texts  dealing  with  the  special 


5.  J.  A.  Stratton,  Electromagnetic  Theory,  McGraw  Hill,  1941, 
section  7.3. 

6.  Abramowitz  and  Stegun,  ed..  Handbook  of  Mathematical  Func¬ 
tions,  AMS  55,  National  Bureau  of  Standards,  1964. 

7.  W.  R.  Smythe,  Static  and  Dynamic  Electricity,  2nd.  ed . , 
McGraw  Hill,  1950. 


8.  Morse  and  Feshbach,  Methods  of  Theoretical  Physics,  McGraw 

Hill,  1953. 


functions  of  mathematical  physics . 6 , 10  For  £  as  a  general  com¬ 
plex  number  not  on  the  cut  (-1  <_  £  <_  1)  the  definition  differs 
somewhat.  The  form  in  equation  8  agrees  with  more  general  defi¬ 
nitions  for  real  m  not  necessarily  an  integer.  In  any  case  we 
only  use  -1  <  £  <  1  in  this  note. 

Similar  to  Stratton  (ref.  5,  section  7.11)  define  two  inde¬ 
pendent  sets  of  vector  wave  functions  which  can  be  used  to  con¬ 
struct  any  divergenceless  vector  field  satisfying  the  vector  wave 
equation  (as  in  equations  6) .  The  first  set  of  vector  wave  func¬ 
tions  is  defined  by 


i5u) 


(n,  m,  ®)  =  V  x  j^r5(£)(n,  m,  ®)ej 


(10) 


where  er  is  the  unit  vector  in  the  r  direction,  and  similarly  for 
other  unit  vectors .  These  vector  wave  functions  have  components 
as 


(n,  m,  ®)  =  0 


M<T(n-  o>  =iIKT9T5T  2‘t)(n'  m' 


U) 


(11) 


V  = 


d  _  U)  /_  e. 

“  (n'  m'  o} 


where  the  last  two  components  can  be  expanded  as 

,in  ' - ,a ' '  -sin  (m<f>) 

cos  (m<f)) 


MqU)  (n,  m,  ®)  =  f^}  (kr)  ~n 


P„ (cos  (0  )  ) 
sin  (0)  m 


MU){n  m  e.  _  fU)(kr)  dPn(COS(8)) 

%  (n'  m'  o]  ”  n  (kr)  “ 


d0 


cos  (m4>) 
sin  (m<j>) 


(12) 


~9~.  Magnus ,  O&erhettinger ,  and  Soni,  Formulas  and  Theorems  for 
the  Special  Functions  of  Mathematical  Physics,  3rd  ed.,  Springer- 
Verlag  New  York,  1966. 

10.  E.  W.  Hobson,  The  Theory  of  Spherical  and  Ellipsoidal  Har¬ 
monics,  Chelsea,  New  York,  1955. 


The  second  set  of  vector  wave  functions  is  defined  by 


NU)(n,  m,  ®)  =  |  V  x  MU)(n,  m,  ®) 


which  has  components 


„U)  _  ex  _  n (n  +  1) 

Nr  (n,  m,  > - jjj— 


-  (n,  m,  ) 


MU),„  „  e.  1  9 2  LrU)<.  m  e.1 

NQ  (n,  m,  Q)  kr  g  (kr )  3 0  [kr“  (n'  m'  o*J 


M  U)  .  _  e.  _ 

Ncj)  ^n'  m'  o  “ 


[krH{£)  (n,  m,  ®)] 


which  can  be  expanded  as 


( ^ )  /  e. 

Nr  (n,  m,  Q) 


f  ^  (kr) 

n  (n  +  1)  - -  P™(cos(9)) 


cos  (met* ) 
sin  (m<j>) 


N0U)  (n'  m'  o> 


.  [krfn 


(kr)  dP™(cos(9))  cos(m4>) 

kr  d0  ~  .  ,  ,« 

sm  (m<p ) 


.  .  [krf ^  (kr>)  P^(cos(0))  -sin(m$)! 

(n'  m'  o>  =  1 - ki - J - sin(6)  m  ,  .. 

Y  cos(m<M 


sin ( 9 , 


where  a  prime  is  used  to  indicate  differentiation  with  respect  to 
the  argument  of  the  spherical  Bessel  function  being  considered. 
Similar  to  equation  13  we  also  have  the  relation 

MU)  (n,  m,  ®)  =  i-  V  x  NU)  (n,  m,  ®)  (16) 


The  N  and  M  functions  are  mutually  orthogonal  on  a  sphere  (of  con¬ 
stant  kr) .  Also  qn  the  unit  sphere  we  have  the  orthogonality  re¬ 
lations  that  two  N  functions  with  different  n  or  m  are  orthogonal, 
or  if  one  is  even  and  the  other  odd  they  are  orthogonal .  The 
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same  results  hold  for  the  M  functions.  For  all  indices  the  same 
we  have 


f  f  3^  (n,  m,  ®)  •  MU)  (n,  m,  ®)  sin(e)ded<J> 

•'A  ~  ^ 


'O  O 


r,  .  o  n  (n  +1)  (n  +  m)  :  \c  U)  /w*l 

=  11  4  W2*  in  +  1  (n  -  ST!  l£n  (kr)J 


f2"  f  3(l>(n,  if  |)  •  fi*8,1  (n,  m,  ®)  sin(e)d0d* 

J  n  •'n 


fo 


r  -i  .  o  ,  n  (n  +  1)  (n  +  m)  ! 

=  Ut  (2n  TIT*  Tfr^5TT 


•  (n  +  l)[f^(kr>]2  +  n[f^[(kr)]2  (18) 

where  5„  „  is  the  Kronecker  delta  function  defined  by 

rnl»m2 

1  for  m^  =  m2 

6  =  US) 

ml'm2  0  for  ^  m2 

With  the  N  and  5  vector  wave  functions  we  can  now  expand  an 
electric  field  with  zero  divergence  as 


co  n 


E  =  EoX)  X  |an(m5  (n'  +  Bn,mN  (n'  m'  o}|  (20) 

n=o  m=o 

where  E0  is  some  constant  with  dimensions  volts/meter  and  where 
an  _  and  Bn  m  are  dimensionless  constants.  Note  that  we  can  also 
sum  over  l  and  over  even  and  odd  functions  but  typically  for  some 
particular  form  of  wave  only  one  l  will  be  needed  and  with  appro¬ 
priate  symmetry  only  one  of  e  and  o  will  be  needed  for  each  type 
of  vector  wave  function.  Comparing  equations  4  and  5  with  the 


relations  between  the  ^  and  N^functions  (equations  13  and  16) 
note  that  we  can  find  H  from  E  by  replacing 

jtU)  /_  _  e>  i  +  ( £ )  ,  e, 

M  (n,  m,  )  -►  ^  N  (n,  m,  ) 


±(SL)  /  ev  i  ;*(£)  /  e. 

N'  (n,  m,  o)  ^  M  ra'  Q) 


giving  an  H  to  go  with  equation  20  as 


E  00  n 

s-^E  Ek, 


n=o  m=o 


mSU)(n,  m,  ®)  +  BnimS(t>(n,  m,  *)|  (22) 


Similarly  if  we  are  given  an  expansion  for  H  we  can  find  E  by 
substituting 

MU)(n,  m,  ®)  +  -iZNU)(n,  m,  ®) 


g(£) (n,  m,  ®)  -*■  -iZ&(£)  (n,  m,  e) 


With  the  N  and  M  functions  we  can  expand  any  electromagnetic 
field  distribution  as  long  as  the  fields  have  no  divergence  which 
requires  p  =  0  in  equations  4.  For  completeness  note  that  the 
wave  equation  for  E  found  from  equations  4  and  5  is 


-V  x  (V  x  g)  +  k22  = 


which  reduces  to 


V2S  -  V(V  •  h  +  k2E  =  0 


As  discussed  by  Stratton  another  set  of  wave  functions  are  needed 
which  we  define  as 

L(£)  (n,  m,  ®)  s  £  VEU)  (n,  m,  ®)  (26) 


which  has  components 
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cos  (m<£ ) 
sin  (mc{> ) 


L^}(n,  m,  ®)  =  f^l)  (kr)P®(cos(0)) 


T  rr,  e\ 

L0  (n'  m'  o) 


f^£)(kr)  dP™ (cos (8 ) ) 


kr 


W 


cos  (m<$>) 
sin  (m<|>) 


(27) 


T  m  e\ 

L4>  (n'  m'  o 


f^  (kr)  P®(cos(0)) 

sin  (0 )  m 


kr 


-si’1  (m<J>) 
cos  (m<j>) 


The  £  and  M  vector  functions  are  related  as 


(n,  m,  ®)  =  krL^  (n,  m,  ®)  *  er  (28) 

The  L  and  M  functions  are  mutually  orthogonal  on  a  sphere  (of 
constant  kr) .  However  the  L  and  N  functions  are  only  orthogonal 
if  at  least  one  of  the  indices  (SL,  n,  m,  ®)  differs.  The  L  func¬ 
tions  are  orthogonal  to  each  other  unless  the  indices  are  all  the 
same,  in  which  case  we  have 


f  f  L^(n,  m,  ®)  •  (n,  m,  ®)  sin  (0)  d0d<)> 


'o  o 


M  .  r  •]  2 7T  (n  +  m)  ! 

m,o  2n  +  1  (n  -  m)  I 


[*n1)<kr>]' 


+  n(n  +  1)1  W  (kr  1 
(kr)2  L  n  J 


(29) 


For  another  reference  concerning  these  spherical  vector  wave 
functions  see  Morse  and  Feshback.^  The  definitions  used  here 
are  similar  to  Morse  and  Feshback  except  that  we  use  elut  v.hereas 
Stratton  and  Morse  and  Feshback  both  use  e”^*-  for  the  harmonic 
time  dependence. 
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III.  Vector  Plane  Waves  in  Spherical  Coordinates 


Now  consider  vector  plane  waves  in  spherical  coordinates. 
Such  a  wave  has  the  general  form 

~y  -y  -i  v  •  r 

F  E  FQue  r  (30) 

where  u  is  some  unit  vector  with  fixed  direction  and  where  the 
propagation  vector  is  given  by 

k  e  (31) 

-y 

where  ei  is  the  direction  of  propagation  of  £he  wave.  Figure  2A 
§hows  a  general  pl^ne^wave  a^  some  position  r  propagating  in  the 
ei  direction  wi^h  E,  H,  and  k  all  mutually  orthogonal  vectors. 

We  assume  that  E  has  a  fixed  polarization  (and  thus  H  also  has  a 
fixed  polarization)  for  purposes  of  illustration.  This  plane 
wave  also  has 

E  —  ZH  x  (32) 


so  that  the  wave  is  propagating  in  the  +ei  direction. 

e^  is  a  fixed  direction_^in  space  determining  the  direction 
of  wave  propagation.  ]£  and  H  are  parallel  to  a  plane  which  is  in 
£urn  perpendicular  to  ei.  Construct  two  orthogonal  unit  vectors, 
e2  and_^e3,  parallel  to  this  plane  as  shown  in  figure  2B.  Now 
since  ei  is  a  fixed  direction  in  space  it  can  be  described  by 
fixed  angles  0i  and  (J>i  in  a  spherical^coordin^te  system  with  re¬ 
spect  to  the  cartesian  unit  vectors  (ex,  ey,  ez)  which  arg  also 
fixed  directions  in  space.  As  shown  in  figure  2§  choose  e2  to  be 
parallel  tg>  the  same  plane  to  \ghich  both  ei  and  eZ->are  parallel; 
also  make*  e2  perpendicular  to  ei.  Finally  choose  e3  perpendicular 
to  both  ei  and  e2,  thereby  making  it  parallel  to  the  x,  y  plane. 

These  new  unit  vectors  are  made  to  form  a  right  handed  system 
so  that  we  have 

-y  -y  -y  ->  -y  -y  ■+  -y  -y 

®1  x  ®2  ~  e3  '  e2  x  e3  =  1  ®3  x  =  e2  (33) 

Not§  that  e2  is  chosen  such  that  for  0  <  0]_  <  tt/ 2  the  polar  angle 
of  e2  is  tt/2  -  8]_.  In  terms  of  the  cartesian  unit  vectors  we  have 
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e,  =  sin  (0, )  cos  (4>t  )  e  +  sin  (0..  )  sin  (<J>,  )  e  +  cos(0,)e 
1  i'  'Tl'  x  1  1  y  1  z 

^  =  -cos  (0^)cos  ( <#> -^ )  ex  -  cos  (0^)  sinf^Je^  +  sinCS^)^ 

-► 

e0  =  sin(<j>,)e  -  cos(4>,)e 
3  1  x  Y1  y 

Cartesian  and  spherical  coordinates  are  related  as 
x  =  r  sin(0)cos  (<J>) 

,y  =  r  sin  (0)  sin  (4>) 

z  =  r  cos  (0 ) 

The  cartesian  and  spherical  unit  vectors  are  related  as 
ex  =  sin  (6)  cos  (4>)er  +  cos  (0 )  cos  (<j> )  e0  -  sin(4>)e^ 

e  =  sin  (0 )  sin  (4>)  e  +  cos  (0 )  sin  (<}>)  eD  +  cos(ij))e. 
y  r  0  <p 

e  =  cos(6)e  -  sin(0)efl 

Z  i  o 

or 

e  =  sin (0 ) cos  (<J>) e  +  sin  (0)  sin  (<J>)  e  +  cos(0)e 

3T  X  y  Z 

eQ  =  cos  (0  )  cos  (4>)  e  +  cos  (0 )  sin  ($)  e  -  sin(0)e 
o  x  y  Z 

=  -sin(<}))ex  +  cos(<j>)ey 

Substitute  for  the  cartesian  unit  vectors  in  equations  .34 
equations  36  and  use  some  trigonometric  identities  to  give 


(34) 


(35) 


(36) 


(37) 


from 
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=  [cos  (61)cos  (6)  +  sin  (0^)  sin  (9 )  cos  (4>  -  4>-,)]er 
+  [-cos  (0^) sin  (0)  +  sin(01)cos  (0)cos  (<f>  -  4>^)]e0 

-  sin  (0^)  sin  (<}>  -  <Jj1  ) 

e2  =  [sin  (G^  cos  (0)  -  cos  (01)  sin  (0)  cos  (<J>  -  4>1)]er 

-  [sin  (6^)  sin  (0)  +  cos  (0  )  cos  ( 0 )  cos  <  4»  -  ?-,)]eQ  (38) 


+  cos  (0^)  sin  ((j)  -  <j)^)e^ 
e3  =  -sin  (0)  sin  (4)  -  4>1)er 
-cos  (9)  sin (4;  -  4>-j_)ee 
-cos  (4)  -  4^)^ 

Note  that  4>  and  4>1  appear  in  the  combination  4>  -  4>i  in  this 
formulation. 

The  general  unit  vector  u  (with  constant  direction)  used  in 
equation  30  can  be  considered  a^  a  linear  combination  of  the  sys¬ 
tem  of  orthogonal  unit  vectors  e]_,  e2 ,  eg.  Now  ei  gives  the  di¬ 
rection  of  propagation  of  the  vector  plane  wave.  For  our  purposes 
we  have  an  electromagnetic  wave  which  is  an  assumed  TEM  plane 
wave.  This  plane  wave  will  be  used  as  the  incident  wave  f$r  cal¬ 
culating  the  sensor  response.  Calling  the^electric  field  Einc 
and  the  magnetic  field  Hinc  then  Einc  and  H^nc  must  both  be  per- 
per.cjicular^to  ej_;  they  can  then  be  formed  from  linear  combinations 
of  e^  and  e2*  In  figure  3  we^illustrate  such  a  plane  wave  inci¬ 
dent  on  a  sphere  centered  on  r  --  0.  The  slot  around  the  sphere 
is  to  be  centered  on  0  =  tt/2  so  as  to  make  the  sensor  symmetric 
witlj  respect  to  the  x,  y  plane.  By  symmetry  if  Einc  were  parallel 
to  e2  and  thus  parallel  to  the  x,  y  plane  it  would  drive  no  net 
current  across  the  gap;  note  we  are  going  to  integrate  the  cur¬ 
rent  all  around  the  gap  to  obtain  the  short-circuit  current  for 
the  sensor.  Thus  we  are  only  interested  in  the  polarization  of 
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the  incident  wave  for  which  E-j_nc  is  parallel  to  and  our  inci¬ 
dent  wave  is  defined  by 


-*■  ->  — ik  •  r 

E.  5  E  e0e 
inc  o  2 


A  o  ■+■  -ik*r 

H .  =  —  e0e 

me  Z  3 


(39) 


Since  one  can  convert  an  electric  field  distribution  to  the  asso¬ 
ciated  magnetic  fiel<|  by  the  transformations  in^equations  21^we 
only  need  to  expand  Einc  in  terms  of  the  M  an£  N  functions;  Hinc 
is  directly  obtainable  from  the  expansion  of  Einc. 


As  a  further  simplification  we  will  later  set  <f>i  =  0  to  make 
k  parallel  to  the  x,  z  plane  as  shown  in  figure  3.  Since  the 
sensor  geometry  is  independent  of  this  represents  no  loss  in 
generality.  Consider  an  r^,  0j_,  <J>i  spherical  coordinate  system 
where  0  jl  and  <J>i  are  the  angles  fixing  the  direction  of  propaga¬ 
tion  of  the  incident  plane  wave  as  above.  The  unit  vectors  for 
such  a  coordinate  system  are 


-* 


e 


r 


1 


e 


1  ' 


(40) 


As  in  Morse  and  Feshback12  we  define  vector  spherical  harmonics 
using  the  spherical  angles  0}./  as 


P(n,  m,  ®) 


(cos  (0^  ) 


cos  (m<{>^) 
sin  (m^) 


B(n,  m,  ®) 


✓ri7n+IT 
(2n+l) sin (0^) 


n-m+1 

"n+T" 


(cos (01) ) 


n+m 

n 


(cos ( 0  ^ ) 


’] 


cos  (m^) 
sin  (mef^) 


.  2  rc(2n+l) 
4^  n (n+1) 


(COS  (01)  ) 


-sin  (m4>^) 
cos  (mcf>^) 


12.  Ref.  8,  Part  II,  pp.  1898,  1899. 


C(n,  m,  ®) 


/n (n+1 ) 
(2n+l) sin (0^ ) 


m(2n+l) 
n (n+1) 


P™  (cos (0, ) ) 
n  1 


-sin  ) 
cos  (m<j>^ ) 


n-m+1 

n+1 


(cos  (6.^ ) 


n+m 

n 


cos  (mtf)^) 
sin  (m<J>^) 


(41) 


The  letters  e  and  o  refer  to  even  and  odd  functions  and  are  asso¬ 
ciated  with  the  upper  and  lower  functions  in  braces  just  as  with 
the  vector  wave  functions  previously  introduced. 

Now  we  introduce  the  result  of  Morse  and  Feshback^  for  a 
dyadic  plane  wave  with  i  switched  to  -i  giving 


(6a.S)e 


-iJc  • 


0 

1 

0 


0\ 

0  e-1^ 


■  E  E  E'2-W'-i)n<2n+1) 

n=o  m=o  e,o 


(n-m) 
(n+m) I 


iP(n,m,®)L(1) (n,m,®) 


£(n,m,®)M^  (n,m,®)  +  iB(n,m,®)N 


/n (n+1) 


(1)  /«  m  e. 
( n  ,  m  ,  ^ ) 


(42) 


The  P,  B,  Q  ve^ctgr  spherical  harmonics  are  functions  of  9f  and  cjjf 
while  the  L,  M,  N  spherical  vector  wave  functions  are  functions 
of  r,  0,  and  4>.  Note  the  use  of  dyadic  notation  where  two  vectors 
are  written  side  by  side  to  form  a  dyadic.  (6a,g)  is  the  identity 
dyadic  which  is  also  written  in  matrix  form  in  equation  42.  A 
vector  plane  wave  is  formed  by  taking  the  dot  product  of  some  con¬ 
stant  vector  on  the  ^eft  of  the  terms  in  equation  42.  Using  a 
general  unit  vector  u  with  fixed  direction  (as  in  equation  30)  we 
have  a  vector  plane  wave  as 


13 .  Ref .  W,  Part  II,  p.  1866. 
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■  £  t,  Ei2-5»,o"-i>”'2«ii 

n=o  m=o  e,o 

[u-C(n,m,®)]M(1)  (n,m,®)  +  i  [u -B  (n,m,®)  ]  N (1)  (n,m,®) 

(43) 

The  problem  of  expanding  a  vector  plane  wave  then  reduces  to 
finding  the  coefficients  expressed  above  as  dot  products. 

Our_^plane  wave^of  interest  (equations  39)  can  be  found  by 
setting  u  gqual  tg  e2  and  then  equal  to  e3  while  using  the  rela¬ 
tions  for  e2  and  e3_>given  in  equations  40  and  setting  4>-l  =  0. 
Considering  first  e2  we  have 

ie2  •  P (n,  m,  ®)  =  -ieQ  *  P(n,  m,  ®)  =  0 


(n-m)  I  .  i: ,  ev,t(l)  ,  e. 
(n+m)  I  iIu*P(n/Ill,0^1'  (n'm'o) 


•  C (n,  m,  ®) 


C (n,  m,  ®) 


m 

n (n+1) 


p^J  (cos  (e1) ) 


0 

1 


(44) 


Pn+l(cos(6l) } 


Thus  only  the  second  set  of  coefficients  for  odd  harmonics  and 
the  third  set  for  even  harmonics  can  bg  non  gero.  Note  that  n  =  0 
is  not  used  because  the  corresponding  M  and  N  functions  are  iden¬ 
tically  zero.  The  last  of  equations  44  can  be  rewritten  using 
identities  for  the  Legendre  functions^  as 


402. 
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.  .  dP  (cos  (9 , ) )  1 

- i -  e„  •  Bln,  m,  ®)  =  .  — 2—™ - - —  (45) 

✓FTimr  2  °  n(n+1)  d0l  0 

The  incident  electric  field  from  equations  39  can  then  be  written 


2.  =  Ee„e 

me  o  2 


-*•  -lk  •  r 


Vlk'r  ■  E  E  k, 

n=l  m=o 


(n,m,o)  +  bn  (n,m,e) 


where 


_w  P™  (cos  (9,  )  ) 

_  r *\  r  i  /  ,;\n+2  2n+ 1  (n  m)  .  n _ 1 

an,m  1  °m,o  n  tn+1)  (n+m) !  sin (0^) 


b  -(25  1 (  i)n+1  2n+1  (n~m> »  dPn(COS(9l)) 

°n,m  u-0m,oJ  l“x;  n (n+1 )  (n+m)!  de1 

Setting  u  =  ^3  we  have  another  set  of  coefficients  as 


ie3  •  £(n,m,®)  = 


P(n,m,o)  =  0 


-  e,  •  £(n,m,®)  =  — — e,  •  J(n,mk) 

✓n (n+1)  /n (n+1)  ^1 


-1  + 


1  1  fn-m+1  „m  , _ ,Q 

iTFry  jH+rnsr-  pn+i(cos(ei)) 


.  n+m  f.  (cos(e  1 

n  n-1  1  J  A 


1 

n  (n+1) 


dP1'1  (cos  (9 , ) )  1 
n  1 
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l  ,  e, 

e3  •  B(n,m,^)  = 


-l 


/n (n+1) 


•  eA  •  B (n,m, _ ) 
/nTn+1)  *1  ° 


(48) 


sin(01)  rTfa+TT  pn(cos  (9X) ) 


The  incident  magnetic  field  from  equations  39  can  then  be  written 
as 


E  .7*-  -► 

H.  =  Z  e"lk*r 
inc  Z  e3e 


(49) 


oo  n 


53e'ik'r  -  E  E  [ian,mSU)(n'm'°>  +  ibn,-S,i'<n 


(1) 


n=l  m=o 


,m,e)] 


We  then  have  the  complete  expansion  of  our  incident  wave. 

Actually,  with  e2e“1K*r  and  e3e“1K*r  expanded  any ^polariza¬ 
tion  of  the  incident  electric  field  (perpendicular  to  k)  can  be 
written  as  a  linear  combination  of  these  two  plane  waves.  The 
associated  incident  magnetic  field  can  be  similarly  written.  Be¬ 
ginning  with  equations  44  we  restricted  =  0.  However,  if  it 
is  desired  to  have  <j>i  /  0  there  are  at  least  two  approaches  to 
generalize  the  present  results.  First,  one  can  simply  leave 
<f>l  ^  0  in  calculating  the  coefficients  as  in  equations  44  and  48; 
this  will  put  terms  like  cos(m<J>i)  and  sin(m<f>i)  in  the  coeffici¬ 
ents.  Or  second,  one  can  take  advantage  of  the  axial  symmetry  of 
the  spherical  coordinate  system  and  consider  <p  -  as  the  azi¬ 
muthal  angle  in  the  present  result^;  in  Jhis  case  one  just  substi¬ 
tutes  <j>  —  in  place  of  <p  in  the  M  and  N  functions  in  equations 
46  and  49.  Which  of  these  methods  if  used  is  a  matter  of  conven¬ 
ience  related  to  the  problem  at  hand.  As  stated  before  we  assume 
the  spherical  sensor  in  the  present  problem  to  have  axial  symmetry 
so  that  we  use  <f>]_  =  0  for  convenience. 

The  expansion  of  e^e-1^’1-  would  be  another  extension  of  the 
present  results  for  cases  where  there  were  source  currents  in  the 
incident  wave  such  that  the  incident  electric  field  had  a  non 
zero  divergence.  ^This,  of  course,  would  be  accomplished  by  sub¬ 
stituting  ei  for  u  and  obtaining  a  set  of  coefficients  analogous 
to  equations  44  or  48. 
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IV.  Short  Circuit  Current 

Having  the  incident  plane  wave,  next  consider  the  short  cir¬ 
cuit  current  from  the  sensor.  For  this  calculation  short  circuit 
the  gap  in  the  sphere  so  that  it  looks  like  a  complete  perfectly 
conducting  sphere  with  no  gap.  Then  we  consider  the  currents  in¬ 
duced  on  the  sphere  by  the  incident  plane  wave.  The  incident  plane 
wave  is  given  by  equations  46  and  49  in  the  previous  section. 

Add  to  this  a  scattered  wave  of  the  form 


Esc  -  EoZ  £  [cn,mS<4)  +  dn,mS(4)(n'"''e>] 

n=l  m=o  J 


°E  £  [ion,mS<4)  (n'm'0)  +  idn,mS<4>  (n'm'e)] 


sc  Z 


n=l  m=o 


The  fourth  kind  of  spherical  bgssel  functions  are  used  to  give  an 
outward  propagating  wave;  the  M  and  N  functions  are  chosen  even 
or  odd  to  match  the  <{>  dependence  of  those  appearing  in  the  expan¬ 
sion  for  the  incident  wave. 

The  boundary  condition  at  r  =  a  is  that  Ee  =  =  0.  From 

the  spherical  vector  wave  functions  in  equations  12  and  15  this 
requires 


a  j  (ka)  +  c  hr  (ka)  =  0 
n,m  n  n,m  n 


[kajn(ka)l •  +  ^  [kahn2) (ka)] ’ 


n,m  ka 
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E  =  E.  +  E  „  , 
inc  sc 


H  =  H.  +  H 

inc  sc 


The  surface  current  density  Js  on  the  perfectly  conducting 
spherical  surface  on  r  =  a  has  0  and  <j>  components  which  can  be 
found  from  the  magnetic  field  just  outside  the  surface  as 


J  =  -H, 

se  *  r=a+ 


oo  n 


■  E  K. 


n=l  m-o 


[kajn(ka)  ]  • 


jn(ka)  Jkah^2)  (ka)j 


h(2)  (ka) 
n 


P“l(cos(G)) 
sin (0) 


m  cos  (m<j>) 


+  ib  '  j„  (ka) 
n,m  n 


[kaj  (ka)]'  (0.  dp  (cos(0)) 

n  u  (2)  n--,\  n 


[kah^2) (ka)j  n 


cos  (m<f>) 


Js,  =  4 

4)  r=a+ 


oo  n 


S  IN. 


n=l  m=o 


[kajn(ka)]'  jn(ka)  ^kah^2)  (ka)j 


h^'  (ka) 


dp“l(cos(0)) 


sin  (m4>) 


[kaj  (ka) ] '  .  P„ 

ib  m  j  (ka)  -  - - A - =p-  h^;  (ka)  — 

n,m  n  f kah  2  (ka)l  n 


P„ (cos  (6) ) 


-m  sin(m<f>) 


2  0 


To  simplify  this  we  use  a  Wronskian  relation^  to  obtain 


(ka) [kajn(ka) ] '  -  jn(ka) 


K2) 


kayn(ka)  [kajn(ka)  ]  '  -  ka jn (ka) [kayn (ka) ] ’ 


(55) 


Using  this  result  in  equations  54  gives 


E  00  n 

rse  ■  ^E  E 

0  n=l  m=o 


n,m  (ka) 2h^2) (ka) 


p”(cos(0)) 

-  -iinirr '  ra 


+  b 


_ -L _ 

n'm  ka£kah^2)  (ka)' 


E  “  « 

Js  ■  srE  E 

<J>  n=l  m=o 


dP™  (cos  (6  ) ) 
x  - 3B - 


dP™ (cos  (0 ) ) 


cos  (m4>) 


(56) 


-a 


n 


n'm  (ka)2h^2) (ka) 


"dTT 


sin  (m4>) 


-  b 


n'm  ka |"kah^)  (ka) 


P^ (cos (0) ) 
sin (0) 


m  sin(m4>) 


With  surface  current  density  on  the  sphere  determined  we  can 
calculate  the  total  current  crossing  a  circle  of  constant  0  on 
the  spherical  surface;  this  is 


15.  Ref .  6 ,  eqn .  10.3.4. 


21 


f  4,  II 

1(0)  =  a  ^infQ)  /  J  d<j> 
JQ  0 


’  2na  .in,.,  bn,o 


dPn (cos (0) ) 


where 


,  „n+l  2n  +  1  dPn<oos(61)) 
bn,o  =  ('1)  n(iT+_rr  - 3?; - 


Rewriting  the  Legendre  functions  with  the  relation 


dP~(cos(0))  , 

~^SS - =  fjlccsfe)) 


we  then  have 


1(0)  =  2 fra  sin (0)  {_i) 


.  x  n+1  2n+l 


kaj^kah^2^  (ka)J 


•  P^  (cos  (0-, ) )  P^  (cos  (0 ) ) 
n  x  n 


For  small  ka  we  have  for  the  spherical  Hankel  functions 


h^2) (ka) 


=  i(2n  -  l)::(ka)"n_1  +  O^(ka)"1^ 


The  double  factorial  function  is  defined  by 


(2n)  : :  =  (2n)  (2n  -  2)  •  *  •  (4)  (2) 


(even) 


(2n  -  1)1!  =  (2n  -  1) (2n  -  3)  • ••  (3) (1) 


(odd) 


with  the  conventions 


O'!  =  1  , 


(-D : : 


1 


(63) 


For  small  |ka|  the  n  =  1  term  in  equation  60  dominates  and  we 
have  as  ka  -*■  0 


1(9) 


3na  ^  (-ika) 


sin(0^) [sin (8) ] 2 


=  -3ira2Eo(a  +  iwe)  sin(0^) [sin(0) ] 2  +  0^(ka)2J  (64) 

Thus  for  low  frequencies  the  short  circuit  current  is  proportional 
to  (o  +  icoe)E0  sin(8i)  which  is  the  z  component  of  the  total  cur¬ 
rent  density.  For  the  case  that  a  =  0  then  the  response  is  pro¬ 
portional  to  iweE0  sin(0i)  which  (in  the  time  domain)  is  just  the 
z  component  of  £).  Setting  0  =  tt/2  so  that  the  position  of  the 
loop  gap  is  on  the  equator  of  the  sphere  we  have  as  ka  -*•  0 


)  sin(01)  +  O^(ka)2^ 
where  the  equivalent  area  of  the  sensor  is 


I  (tt/2)  =  -AeqEQ(a  +  iwe] 


(65) 


A  »  3ira" 
eq 


or  as  a  vector 


-  3Tra2e 
eq  z 


so  that  we  can  write  as  ka  0 
I(tt/2)  =  -(a  +  iwe)E 


me 


,  •  %  + 
?=o  eq 


0  ^(ka)  2^ 


(66) 


(67) 


(68) 


For  convenience  we  define  a  short  circuit  current  transfer 
function  as 


T(81)  =  I(tt/2)  [~AeqE0(c  +  iwe)]"1  (69) 

so  that  as  ka  0  we  have  T  -*■  sin(8]J  which  is  the  ideal  low- 
frequency  angular  dependence  of  the  sensor.  Then  we  have 
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T<9l>  =  l?i('i>n  ^  (ka)2^(2,(ka)J'  Pn(C°S,e1»Pn(0'  <7°> 

Some  special  values  for  the  Legendre  functions  are1^ 


P>) 


0  for  n  +  m  odd 
n+m 


(71) 


(-1) 


2  (n  +  m  -  1)  !  I 
(n  -  m)  !  ! 


for  n+m  even 


giving 


°°,2 


T(01>  =  E  Jn 


2  2n+l 


n:  ! 


T  n(n+1)  (n-1)  : !  (ka)2  kah^2*  (ka) 

_ 


rPjfcoste^)  (72) 


where  the  second  index  above  the  summation  sign  is  the  increment 
which  is  added  to  the  value  of  n  to  obtain  the  next  value  of  n 
for  the  summation.  Here  only  odd  n  result.  For  convenience 
define 


W 


n 


0)  =  (~i)  t 


2 

3 


2n+l 
n (n+1) 


(ka) 


"kah^2)  (ka)  ’ 


pJ-(cos(6  ))pJ-(0)' 
n  in 


(73) 


For  0  =  tt/2  and  odd  n  this  is 


W 


tt/2)  = 


2n+l 
n (n+1) 


n : 


(n-1) 


(ka) 


kah(2)  (ka)l  ' 

L  n  J 


P  (cos  (0. ) ) 
n  i 


(74) 


For  our  present  case  of  interest  just  call  this  so  that  the 
short  circuit  current  transfer  function  can  be  written  as 


T(0X) 


n=l 


(75) 


Note  that  as  ka  ^  0  we  have  T  -*•  Tj.  +  sin(0]L)  . 


16.  Ref.  7,  p.  153. 
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The  short  circuit  current  transfer  function  is  plotted  in 
figure  4  as  a  function  of  ka  with  a  =  0.  Note  that  we  plot 
T/sin(0jJ  for  several  values  of  0^.  This  is  so  that  all  the  curves 
tend  to  one  for  low  frequency  and  we  can  observe  at  what  frequency 
(i.e.  what  ka)  they  start  to  spread  from  one  another.  For  compar¬ 
ison  Ti  from  equation  74  is  included.  T]_  will  be  used  later  in 
calculating  some  of  the  sensor  response  functions.  In  figure  4 
the  magnitude  and  phase  of  T/sin(8^)  are  plotted.  For  convenience 
the  phase  is  plotted  as  arg(T)  -  ka  which  corresponds  to  multiply¬ 
ing  T  by  e“i^a.  This  is  the  same  type  of  display  as  used  in  ref¬ 
erence  2.  Note  that  the  magnitude  of  T/sin(0]_)  peaks  up  slightly 
at  around  ka  =  .7.  Also,  between  ka  =  1  and  ka  =  2  the  response 
functions  for  various  values  of  0^  begin  to  spread  apart;  some  ka 
roughly  between  1  and  2  can  then  be  considered  the  upper  frequency 
for  which  the  sensor  maintains  a  response  proportional  to  sin  (02.) 
which  can  be  considered  the  ideal  angular  dependence  of  the  sensor. 
This  range  of  ka  is  also  where  the  frequency  response  (for  the 
short  circuit  current)  starts  to  fall  off  for  higher  frequencies. 

T^/sin(0^)  is  a  useful  function  to  use  for  some  of  the  later 
response  function  calculations  because  it  is  the  first  term  in  the 
expansion  of  T/sin(0i)  and  because  it  has  no  dependence  on  0].. 

The  Hankel  function  for  n  =  1  is  just 


(76) 


Then  T^  is  just 


sin  (0^ 


_ -i _  _  _ 

(ka)2  kah^2^  (ka)  1  +  ika  -  (ka)2 


ika 


(77) 


and  we  also  have 


_  -ika 
T1e 

sin (@1) 


-1 


=  [1  +  ika  -  (ka)  ] 


(78) 


T1/sin(0i)  has  a  peak  magnitude  of  2//T  =  1.155  occurring  at 
ka  =  1//7  -  .707. 
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V.  Admittances 


Nov;  consider  the  admittances  when  the  sensor  is  driven  at  the 
gap.  As  illustrated  in  figure  5  the  loop  gap  is  centered  on 
0  =  tt/2  and  has  an  angular  width  of  2\pQ.  There  is  a  voltage  Vgap 
uniformly  distributed  around  the  gap.  Associated  with  Vgap^there 
are  three  surface  current  densities  which  are  parallel  to  e$. 

These  are  Jsext'  Jsint'  and  Jsc  an<^  are  associated  respectively 
with  fields  external  to  the  sphere,  fields  internal  to  the  sphere, 
and  currents  into  the  cables  or  other  transmission  lines  loading 
the  gap.  Taking  the  conventions  for  the  directions  of  these  cur¬ 
rents  as  indicated  in  figure  5  we  have  three  admittances  to  define, 
namely 


J  J 

v  _  0  _  sint  v  0  _  sint  v  0  „  c 

Yint  ~  2lTa  V  ___  '  Yext  =  2lT  V_.__  '  Yc  “  2ira  V 


(79) 


gap 


In  normalized  form  we  also  define 


gap 


gap 


Yint  ZYint  '  Y 


ext 


ZY 


ext  ' 


y  E  ZY 


(80) 


For  the  numerical  calculations  we  take  a  =  0,  as  before.  Also  de¬ 
fine  a  normalized  cable  conductance  as 


r 


c 


(81) 


where  Zc  is  the  net  cable  impedance  (resistive)  loading  the  gap. 
Since  we  have  a  =  0  for  the  numerical  calculations  then  rc  >_  0  for 
such  calculations  and  rc  is  a  constant  which  we  can  specify  para¬ 
metrically. 


A.  Boundary  Conditions  at  Gap  in  Sphere 

Here  we  use  approximately  the  same  quasi  static  approxi¬ 
mation  for  the  electric  field  distribution  in  the  gap  as  was  used 
in  references  1  and  2.  Defining 

4)  s  j  -  0  (82) 


we  have 
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r=a 


(83) 


E 


0 


with 


fE  (?) 


1  ?  -1/2 

i[l  -  zn  for  U|  <  1 

0  for  |  ?  |  ^  1 


(84) 


This  field  distribution  has  the  proper  form  of  singularity  at 
the  edges  of  the  gap,  considering  the  assumption  of  a  perfectly 
conducting  spherical  shell  of  zero  thickness.  The  presence  of 
cable  connections  across  the  gap  will  of  course  distort  this  field 
somewhat  near  the  connections.  However,  of  the  various  choices 
for  the  gap  field  one  might  choose  the  above  choice  seems  to  be  a 
reasonable  one.  Note  that  we  assume  ipQ  <<  1  for  these  calculations. 
Later  we  change  to  another  field  distribution  f£  which  closely  ap¬ 
proximates  f £  for  small  ip.  This  other  distribution  will  be  more 
convenient  for  the  evaluation  of  a  certain  integral. 


B.  Internal  Admittance 

To  calculate  the  internal  admittance  we  expand  the  fields 
inside  the  sphere  in  terms  of  the  spherical  vector  wave  functions. 
The  boundary  conditions  at  r  =  a  are  taken  as  independent  of  (j>  so 
that  m  =  0  for  the  wave  functions.  There  are  three+non  zero  field 
components:  Er,  E0,  H<{> .  The  fields  are  finite  at  r  =  0  so  that 

the  spherical  Bessel  functions  of  the  first  kind  are  used.  The 
fields  are  then  expanded  as 


®int  =  e> 

n=o 


°-  e» 

n=o 


(85) 


where  E^  is  some  constant  with  dimensions  volts/meter.  The  non 
zero  components  of  the  spherical  vector  wave  functions  for  m  =  0 
from  equations  12  and  15  are 
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{£)  {£)  dP  (cos(0))  /n\  , 

Mj  '(n,0,e)  =  -f^;  (kr)  -■  ^ (kr)  (cos  (0) ) 


f  ^  (kr) 

NT;  (n,0,e)  =  n{n  +  1)  -2-^ -  P°(cos(0)) 


N^£) (n,0,e)  = 


krf  (Z)  (kr) 


]’  dPn 


(cos  ( 0 ) )  krfjil)  (kr) 


]'.i 


Pn(COS  (0)  ) 


The  boundary  conditions  for  Eq  on  r  =  a  are  given  by  equa¬ 
tions  83  and  84.  Writing  out  Eg  on  r  =  a  gives 

"  [kaj  <ka)T  L 

Ee  =  El£“n - h - Pn(cos<9>>  <87) 

r=a  n=o 

Multiplying  both  sides  by  P„(cos(0))  sin (6),  integrating  over  0 
from  0  to  it,  and  using  the  orthogonality  of  the  Legendre  functions 
gives^7 

<88> 

o  o  o 

For  the  integral  we  define  a  convenient  term 

AA  =  r  rfEfrVi(oosl0))  sin<e>d« 

ro  o  \  o/ 


=  i-  f  fi(S-)pi(sin(*))  cos  (^)di^ 

*o  */-tt/2  E' V  n 


=  t1  '  ^  J  Pn(Sin<*))[1  '  Si 


sin  (^)] 


17.  Ref.  6,  eqn.  8.14.13. 
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Since  we  assume  \pQ  <<  1,  then  for  Ji{)|  £  \pQ  we  can  use  sin  (4))  =  ip 

to  write 


i  Ar  /*  \2i  *"  2 1/2 1 

AA  *  ^  U  L1  '  J  11  '  1  P»('t>d*  =  A" 


-1/2 


Then  let  z;  =  'P/^Po  to  give 
.1 


if  2  T  2  1 

Aa-?j[1-5  l  L1  "  ft*] 


Since  the  form  of  fE  is  only  approximate  we  take  An  instead  of  An 
to  calculate  the  coefficients.  This  integral  is  treated  in  ap¬ 
pendix  A.  Using  the  approximation  good  for  small  ipQ  we  then  have 
for  our  coefficients 


gap  2n  +  1  ka  , 

“n  ~  aE^~  2nJn~+~TT  TkaJ^lklJr  An  (92) 

In  obtaining  the  result  of  equation  91  for  An  the  approxima¬ 
tion  sin  (40  =  4>  was  used  for  small  ip.  This  is  equivalent  to  alter¬ 
ing  slightly  the  field  distribution  fE  to  one  which  we  call  fE* 

To  see  this  let  ipo?  =  sin  (40  in  equation  91  giving 


A_  = 


/•arcsin  (4^_ ) 

/  \l 

.  /  sin  (40  \\ 

'-arcsin  (ip  )  L 

1  (  *o  >  J 

-1/2 


cos  (40  (sin  (40 )  cos  (40  d4> 


1  /*1T/ 2  1 

=  -r- /  f  p";  (sin(40  )  cos(4>)d4» 
VnJ-i r/2  E  11 


^0^/2  “  11 
where  we  have  defined 


,“1/2 


■  r  •  2 1 

i  1  -  /S1P  w.)  cos  (40  for  |4>|  <arcsin(4)Q) 

L  J 


0  for  1 4>  |  arcsin(4)o) 
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The  voltage  across  the  gap  is  given  from  equation  83  (except  using 
f^)  by  the  integral 


V. 


V  /-it/2 
-2*2/ 

a*o  *'-tt/2 


fEad^ 


(95) 


Substituting  from  equation  94  and  letting  =  sin(\p)  we  have 


V 


o 


V 

gap 


(96) 


Thus  this  field  distribution  is  still  exactly  consistent  with  Vgap 
as  the  gap  voltage.  The  angular  gap  half  width  is  arcsin(ij>0)  for 
this  distribution  and  is  very  nearly  \pQ  for  small  t fi.  We  thus  take 
fE  as  the  field  distribution  across  the  gap  instead  of  the  simpler 
form  given  by  f£- 

The  surface  current  density  associated  with  the  internal  ad¬ 
mittance  is  taken  at  0  =  ir/2  for  convenience  so  that  we  have 


int 


-Hintj 


r=a 

0=ir/2 


(97) 


From  equations  85  we  then  have 


Sint 


(98) 


n=o 


From  equation  71  we  then  have 


int 


- 


n+1 

2 


n=l 


fri" "-'T'rrr  anjn(ka) 


(99) 


Note  that  only  odd  n  contribute  to  this  sum. 


From  equations  79  and  80  the  normalized  internal  admittance 
is  then 
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(100) 


y.  =  iiaZ 
yint  V  s 


„  2TraE,  ®f2 

—  j  =  - - L  Y'ip1(0)a  j  (ka) 

s.  .  V  L-j  n  nJn 

gap  xnt  gap  n=l 


Substituting  for  an  gives 


*int  -  )L>  HTE~7TT  pnIQ)An 


kajn(ka) 


(101) 


From  appendix  A,  equation  A21,  we  have 


An  *  pi<°>p(-  *TT’  l!  1;  *o) 


(102) 


where  F  is  a  hypergeometric  function  given  by 


^  (-  ¥)  (!) 


n+1  n  ,  .  2\  v-' 

F(-  -j-,  j;  1;  *1  =  X, 

'  '  rt—r 


(qD2 


and  where  (a)q  is  a  Pochhammer  symbol18  defined  by 


(a)Q  =  1 


(a)^  =  a(a+l)(a+2)  •••  (a+q-1)  for  q  =  1,  2, 


(103) 


(104) 


Note  that  the  series  in  equation  103  has  only  a  finite  number  of 
terms  giving  a  polynomial  function  of  ij>0.  This  is  convenient  for 
the  numerical  calculations.  The  normalized  internal  admittance 
can  then  be  written  as 

„  -  V-2,-  2n+l  r„l,„,l2„/  n+1  n.  ,2)  kaVka) 

pint  If  n  (n+1)  lPn^J  p\  2  ’  2'  '  /  [kaj  (ka)  J  ’ 


_  V'  4 _  2n+l  T  nl!  r„/_  n+1  n  ,  ,2\ 

2^17r  n (n+1)  L(n-i)  1 ! J  F  \  2  '  2;  1;  ^o  ) 

n=l 


ka jn (ka) 


[kajn (ka) ] 7 


(105) 


"I  Ref .  6 ,  eqn .  6.1.22. 
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The  asymptotic  properties  of  this  series  for  large  n  are  treated 
in  appendix  B. 

The  normalized  internal  admittance  is  plotted  in  figure  6  as 
a  function  of  ka  with  a  =  0  for  several  values  of  ip0.  Note  that 
Yint  has  only  an  imaginary  part.  For  convenience  yint/ha  is 
plotted  in  figure  6;  yint/ka  is  an  imaginary  function  for  real  ka 
which  we  are  using.  As  ka  0  yint/ka  tends  to  a  constant;  the 
numerically  determined  coefficients  are  listed  in  table  1. 


Table  1. 


4* 

o 

yint 

ika 

001 

14.22 

.01 

9.62 

.1 

5.10 

Asymptotic  form  of  Y^nt  for  small  ka 


As  one  would  expect  for  small  ka  this  admittance  represents  a 
capacitance.  As  ka  is  increased  y^nt  has  singularities  correspond¬ 
ing  to  the  zeros  of  [kajn (ka) ] ' .  The  first  two  singularities  oc¬ 
cur  at  ka  -  2.744  and  ka  =  4.973. 


C .  External  Admittance 

The  calculation  of  the  external  admittance  follows  the 
same  development  as  the  internal  admittance.  The  fields  outside 
the  sphere  are  expanded  as  outward  propagating  waves  in  the  forms 


ElZXS(4)  (n'°'e) 

n=o 


fiEi6nS(4)  (n'°'e) 

n=o 


(106) 


The  only  difference  on  the  boundary  r  =  a  between  these  fields  and 
those  of  equations  85  is  the  replacement  of  jn  (ka)  by  hj\2)  (ka)  . 
Analogous  to  equation  92  we  then  have  the  expansion  coefficients 
given  by 


gap  2n  +  1 


'n  aEx  2n(n  +  1)  [kah 


ka 

ttt 


n 


(ka)] 


i  A 
1  n 


where  only  odd  n  are  of  interest, 


(107) 
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(110) 

This  series  is  very  similar  to  the  one  in  equation  105  and  its 
asymptotic  properties  for  large  n  are  treated  in  appendix  B. 

The  normalized  external  admittance  is  plotted  in  figure  7  as 
a  function  of  ka  with  c  «  0  for  several  values  of  ^0.  In  this 
figure  the  real  and  imaginary  parts  of  yext  are  plotted,  again  in 
the  form  yext/ka  which  tends  to  a  constant  for  small  ka;  the  nu¬ 
merically  determined  coefficients  are  listed  in  table  2. 
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4* 

ro 

v 

J  ext 

ika 

.001 

17.36 

.01 

12.75 

.1 

8.04 

Table  2.  Asymptotic  form  of  Yext  for  small  ka 

Yext  is  also  a  capacitance  for  small  ka.  Note  that  as  ka  is  in¬ 
creased  Yext  does  not  have  singularities  as  was  the  case  with  yint* 
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VI .  Frequency  Response  Characteristics 


Now  that  we  have  the  short  circuit  current  and  the  admittances 
in  suitably  normalized  forms  we  combine  these  results  to  calculate 
the  frequency  response  characteristics.  First  define  a  response 
function  including  only  the  admittances  as 


R 


yint  +  yext 


TT  ■  11  +  rc<yint  +  W1 


-1 


(111) 


This  is  plotted  as  a  function  of  ka  in  figures  8  and  9  for  two 
values  of  i|;0  (.01  and  .1)  with  c  =  0.  For  each  graph  several  val¬ 
ues  of  rc  are  used.  Note  for  rc  >  0  that  Ry  has  zeros  at  the  sin¬ 
gularities  of  yint*  As  one  would  expect  decreasing  rc  maintains 
Ry  as  a  flat  response  characteristic  out  to  larger  values  of  ka 
which  represents  higher  frequency  response. 

Including  the  short  circuit  current  transfer  function  from 
equation  75  we  have  the  response  function 


T(0  )  TO.) 

R^91^  ~  sin (0.  )  Ry  ~  sin  (0.  )  ^  +  rc  ^yint  +  yext^ 


(112) 


Another  convenient,  but  somewhat  artificial  response  function  is 
given  by 


T1  T1 
R1  ~  sin (0. )  Ry  ~  sin  (0. ) 


[1+r  (y.  .  +  y  .)] 
c  1  mt  ■'ext 


-1 


(113) 


This  last  response  function  uses  only  the  first  term  in  the  expan¬ 
sion  for  T(0i)  and  is  independent  of  @i. 

In  figures  10  and  11  we  have  Ri  plotted  as  a  function  of  ka 
for  two  values  of  4>0  (.01  and  .1)  and  several  values  of  rc.  Note 
that  the  factor  Ti/sin(0i)  peaks  up  the  frequency  response  in  the 
vicinity  of  ka  =  1,  partially  compensating  for  the  rolloff  with 
frequency  associated  with  Ry. 

Based  on  R]_  we  define  an  upper  frequency  response  as  the  min¬ 
imum  value  of  ka  for  which 

|R,  I  =  —  (114) 

/2 
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This  value  of  ka  is  plotted  as  a  function  of  rc  in  figure  12  for 
two  values  of  tpQ.  Frequency  response  is  increased  by  increasing 
ipo  and  by  decreasing  rc.  It  may  not  be  desirable  to  decrease  rc 
too  far  because  as  rc  +  0  no  power  is  delivered  to  the  output. 

In  figures  13  and  14  we  have  plotted  R  as  a  function  of  ka 
with  a=0  for  various  values  of  Oq  with  i|>0  =  .1  for  two  specific 
values  of  rc.  These  two  values  are  rc  -  .1327  and  rc  -  .2654 
which  correspond  to  Zc  =  50  ft  and  Zc  =  100  ft  respectively  if  the 
media  inside  and  outside  the  sphere  are  assumed  to  have  the  same 
constitutive  parameters  as  free  space  so  that  the  wave  impedance 
is 

Z  =  ZQ  -  376.7  ft  (115) 

Using  Ri  as  an  average  of  R  to  remove  the  dependence  on  Og,  the 
frequency  responses  for  these  two  cases  as  defined  by  equation  114 
are  given  by  ka  =  .60  and  ka  -  .30  for  Zc  =  50  ft  and  Zc  =  100  ft 
respectively. 
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VII .  Summary 


In  this  note  we  have  developed  frequency  and  angular  response 
curves  for  the  hollow  spherical  dipole  with  a  uniformly  resistively 
loaded  equatorial  slot.  For  these  calculations  the  media  inside 
and  outside  the  sphere  were  assumed  to  have  the  same  permittivity 
and  permeability  and  zero  conductivity.  For  low  frequencies  the 
response  of  this  kind  of  sensor  is  proportional  to  one  component 
of  the  displacement  current  density  (which  is  also  called  the  time 
rate  of  change  of  the  displacement  vector) . 

As  an  extension  of  the  present  calculations  one  might  consider 
the  case  that  the  external  medium  was  a  linear  conducting  medium. 

In  this  case  if  the  resistance  due  to  the  signal  cable  loading  (Zc) 
were  small  compared  to  the  resistance  loading  due  to  the  conduct¬ 
ing  medium,  then  for  low  frequencies  this  sensor  would  have  a  re¬ 
sponse  proportional  to  the  total  current  density  (conduction  plus 
displacement)  in  the  external  medium. 

Perhaps  the  design  considered  in  this  note  can  be  extended  to 
the  case  of  several  slots  around  the  sphere  at  different  values  of 
0 .  This  would  be  a  multi-gap  spherical  dipole  analogous  to  the 
multi-gap  cylindrical  loop.  Another  extension  of  the  present  cal¬ 
culations  would  be  to  loops  with  spherical  geometries.  Perhaps 
some  of  these  and  related  topics  can  be  considered  in  future  notes. 
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A.  PLANE  WAVE  WITH  FIXED  POLARIZATION 


B.  UNIT  VECTORS  FOR  PLANE  WAVES 


FIGURE  2.  VECTOR  PLANE  WAVES  IN  SPHERICAL  COORDINATES 
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Hjnc  IS  POINTING 
OUT  OF  THE  PAGE. 

y  IS  POINTING 
INTO  THE  PAGE. 


FIGURE  3.  PLANE  WAVE  INCIDENT  ON  SPHERE 
CROSS-SECTION  VIEW  WITH  ^=0 
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2  0,/ir  =.25 


sin(0j) 


2  x  10' 


A.  MAGNITUDE  OF  T/sin(0.)  WITH  0,  AS  A  PARAMETER 


arg(T)-ka 


arg(T|  )-Ka 


20,/ tt  =.25 


B.  PHASE  OF  T  WITH  0,  AS  A  PARAMETER 


FIGURE  4.  SHORT  CIRCUIT  CURRENT  TRANSFER  FUNCTION 


i 


FIGURE  5.  CURRENTS  ON  SPHERE  FOR  ADMITTANCES: 
CROSS  SECTION  VIEW 
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B-  ^m[yext]^Q 


FIGURE  7.  NORMA 


FIGURE  8.  EFFECT  OF  ADMITTANCES  ON  RESPONSE:  *  =01 

•  0 
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ka=4.973 


0 
-I 
-2 

arg(R|)-ka 

-3 
-4 
-5 

I0~2  I0'1  ka  10°  10 1 

B.  PHASE  OF  R,  WITH  rc  AS  A  PARAMETER 


FIGURE  10.  ARTIFICIAL  RESPONSE  CHARACTERISTS:  t|/o  =  .0! 
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ko  =4.973 


I0'2  I0'1  re  10°  I01 


FIGURE  12.  DEPENDENCE  OF  FREQUENCY  RESPONSE  ON  CABLE  IMPEDANCE 
WITH  t0  AS  A  PARAMETER:  |R,|=  ufi 
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I0"2  I0"1 

B.  PHASE  OF  R  WITH  0 


FIGURE  13.  RESPONSE  CHARACTERI 


A.  MAGNITUDE  OF  R  WITH  0,  AS  A  PARAMETER 


I0~2  IO"1  ka  10°  10 1 

B.  PHASE  OF  R  WITH  0,  AS  A  PARAMETER 


FIGURE  14.  RESPONSE  CHARACTERISTICS:  ^=.1,  rc  =  lOOfl  /  Z0  -.2654 
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Appendix  A:  Calculation  of  An 


In  equation  91  we  have  a  function  which  is  expressed  as  an 
integral  given  by 


A 


n 


ty0C>dc 


(Al) 


This  function  appears  in  the  coefficients  of  the  field  expansions 
used  to  calculate  the  internal  and  external  admittances  of  the 
sensor.  The  purpose  of  this  appendix  is  to  calculate  this  integral. 


Define 


A-l 


[i  -  ri 


-i 


Rearrange  this  as 


}  -  <&2] 


-iM/2 


Pv<V)d? 


f  2  X-1r  2  21^^  u 

“1  -  J  [1  -  c2]  [l  -  4^2]  P>0C)dC 

-1 


+  fn-  ?2]x_1[i  - 


■ 


A-l 


ti  -  ri 


t1  -  &] 


M/2 


/ 


A-l 


Pu(-^  c)d? 
v  To 


[l  -  *o'2]',/2<<V>dt 


+  [i  -  ri 


giving 


(A2 ) 


(A3) 
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1A.  Ref.  6,  eqn.  15.4.23. 
2A.  Ref.  9,  p.  53. 
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1/2-m+I 


^(1 


■  TLLZ21  _  F/_  H±v  l-P+v.  1.  42  2\ 

-  \  -  M  -  ^  + 1)  (  2'  2  '2  °  j 


(A8) 


This  then  gives 


rA 

y  V  \T 

.  /,  y  ,  v  \ 

Ml  ' 

2  ~  2/ 

l1  "  2  +  2/ 

[1.?2jX'1f/_  »+v  1-M+v 


(- 


The  hypergeometric  function  is  given  by 

"  (a)  (8)  q 

F(a,  6;  y;  z)  =  F(B,  a;  y;  z)  =  ]>^  — fyT — ^  qT 

q=0  9 

where  the  Pochhammer  symbol^A  is  defined  by 


(a)o  =  1 


(a)  =  a(a  +  1) (a  +  2)  •••  (a  +  q  -  1) 

and  provided  a  is  not  a  negative  integer  or  zero 

(a)  =  Ha  +  q) 

lcUq  Ha). 

Then  we  have  the  integral 

h2  HJV. 


•  b 

(A9) 


(A10) 


(All) 


(A12 ) 


E 

q=0 


(-  ♦?«  ,1  .  X-l 


_ C - 

(I) 


a 

q1  ^ 


[i-c2] 


C2qd^ 


(A13) 


3A.  Ref.  6,  eqn.  6.1.22, 
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This  leaves  us  with  the  integral  (in  which  we  let  £  =  as 
rl  A-l 

H3  =  /  [1  -  ri  c  q<u 

•'o 


where  B  is  the  well  known  beta  function  and  we  need  Re [A]  >  0. 
Recombining  these  terms  gives 


(A16 ) 


This  result  for  Hi  is  very  similar  to  a  formula  in  a  standard 
reference.^  However,  in  the  course  of  our  present  investigations 
we  found  this  formula  to  be  in  error,  most  likely  through  a  mis¬ 
print.  This  led  us  to  do  the  above  derivation  to  find  the  correct 
formula.  We  would  like  to  thank  Prof.  Fritz  Oberhettinger  of 
Oregon  State  University  for  our  discussions  with  him  on  this 
problem. 

4A.  A.  Erdelyi,  ed.,  Tables  of  Integral  Transforms,  Vol.  2, 

McGraw  Hill,  1954,  p.  318,  eqn.  31. 
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In  equations  A2  and  A16  let 


y  =  l  , 


v  =  n 


Then  we  have 


n  7T  “1 


(A17) 


2rfi 


W  F(.  2+i  n.  2\ 

:$\rfei\FV  r-r1'  *o) 


(A18) 


Since  (n+l)/2  is  a  positive  integer  we  have 


r/Wl\  _  /n-l\,  _  (n-1)  I ! 

\ 2 )~  \ 2  r  '  aEi 


r  -  ^ 


-(-!)(-  5* 


(-D 
.  2 


(A19) 


This  gives 


An  =  (-1) 


nil  _  /  n+1  n  ,  .  2\ 

nm  F\;  ~r •  v  1;  *0) 


Using  equation  71  this  can  be  written 
*»  ■  p>)F(-  1;  *0) 


Writing  out  the  series  we  have 


(A20) 


(A21) 
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(A22 ) 


Note  that  since  (n+l)/2  is  a  positive  integer  the  series  in  equa¬ 
tion  A22  has  only  a  finite  number  of  terms,  making  An  a  polynomial 
in  'f'o. 
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Appendix  B:  Properties  of  Admittance  Series  for  Large  n 


The  most  difficult  parts  of  the  sensor  response  functions  to 
calculate  are  the  internal  and  external  admittances  because  of  the 
slow  convergence  of  the  series.  In  this  appendix  we  consider  the 
behavior  of  these  series  for  large  n. 


From  equation  105  the  normalized  internal  admittance  is 


'  1 2 


int 


-  E 


ITT 


n=l 


2n+l 
n (n+1) 


n+1 
2  ' 


x  kajn(ka) 

)  [kaj~ (ka) J  ' 


(Bl) 


and  from  equation  110  the  normalized  external  admittance  is 


■^ext 


00  /  2 

E  (-i),r 

n=l 


2n+l 
n (n+1) 


n+1 
2  ' 


(B2) 


where  from  equation  71  we  have 

n+1 

P^(0)  =  (-1)  2  (n^j,  ,  for  n  odd  (B3) 

Define  the  individual  terms  in  the  summation  for  yint  as  yintn 
and  in  the  summation  for  yext  as  yextn  so  that  equations  Bl  and  B2 
can  be  rewritten  as 


“> ,  2 

^int  -  ^int 
n=l  n 


and 


“,2 

^ext  ~  S  ^ext 
n=l  n 


Define 


n-^  H  n2  -  2  >  0 


(B4 ) 


(B5) 


(B6 ) 


where  ni  and  n2  are  odd  integers.  Then  write  the  normalized  ad 
mittances  as 
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In  this  appendix  we  consider  the  individual  terms  in  the  sums  and 
the  error  terms  for  large  n. 

Consider  the  parts  of  the  individual  terms  in  the  admittance 
sums.  Define  an  integer  as 

N  =  -  t— -  (B9) 


where  only  odd  n  >  0  are  of  interest.  For  the  hypergeometric  func¬ 
tion  we  have 


q=o 
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(BIO) 


N  (-N)(N  -  4) 

T  — 

^  (q  DZ 


Now  this  particular  hypergeometric  function  can  be  written  as  a 
Jacobi  polynomial  usinglB 


F(-N,  a+l+B+N;  a+1;  x)  =  P^a'g)  (1  -  2x) 


which  gives 


(°'~  t)/  s'! 

Fn  =  V  /{1  -  2*o) 


The  Jacobi  polynomials  can  be  written  as^B 


P^a'e) (x)  =  2 


ina 


(x  -  1)  q(x  +  l)1 


where  binomial  coefficients  are  given  by 


(Bll) 


(B12) 


(B13 ) 


(B14 ) 


/N+a\  _  r  (N  +  a  4-  1)  =  (N  4  a  -  q  4  1), 

\  q  )  q1.  T(N  +  a  -  q  +  1)  q! 

/N+6\  _  r(N  f  6  +  i)  _  (q  +  g  +  1)n- 
\N_q/  (N  -  q)  :r  (q  +  6  +  1)  (N  -  q)  : 


The  hypergeometric  function  of  interest  can  then  also  be  written 


Fn  =  2 


TbT  Ref .  6 ,  eqn .  15.4.6. 
2B.  Ref.  9,  p.  211. 
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■  J)(*r  <■  -  *!)• 


(B15) 


Note  that  the  Jacobi  polynomial  in  equation  B13  is  a  polynomial 
function  of  both  a  and  0  as  well  as  its  argument;  the  same  is  true 
of  the  hypergeometric  function  in  equation  Bll.  The  range  of  a 
and  0  are  often  restricted  to  a  >  -1  and  3  >  -1  in  order  to  make 
the  weight  functions  integrable . 3B  However  we  are  not  concerned 
here  with  the  orthogonality  properties  of  the  Jacobi  polynomials. 
Since  the  Jacobi  polynomials  are  polynomial  functions  of  a  and  £, 
as  is  the  hypergeometric  function,  then  both  series  are  identical 
for  all  a  and  0.  Thus  B  =  -{3/2)  is  allowed  for  the  series  repre¬ 
sented  in  equation  B15. 

With  the  hypergeometric  function  represented  as  a  Jacobi  poly¬ 
nomial  we  look  at  the  asymptotic  form  for  large  N  (implying  large 
n)  of  the  Jacobi  polynomials.  As  N  »  for  0  <  5  <  tt  we  have^B 


P^a'8)  (cos  (O  ) 


cos  (n  +  i(a+3+1))s  -  J(l+2a)  ( 

- - 2 - L - 4 - _j.  +  0\n  2) 


(big) 


cos  (c)  =  1  -  2t|/ 


(B17) 


giving 


in(§)  - 


-  cos  (?) 


(§)  -  p 1  r4j-]1/2  - b  -  *o]: 


(bib) 


3b.  A.  Erdelyi,  ed..  Higher  Transcendental  Functions,  Vol.  2, 
McGraw  Hill,  1953,  p.  168. 

4B.  Ref.  9,  p.  216. 
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Applying  this  to  equation  B12  for  0  <  <  1  gives  as  n 


-*•  GO 


(B19) 


Thus  as  n  ■+  «  the  envelope  of  Fn  falls  off  proportional  to  n“(l/2). 
Also  note  as  4>0  -*•  0  we  have 

arccos^l  -  2^  =  2\pQ  +  (B20) 

As  a  check  Fn  was  calculated  from  equation  BIO  and  the  asymptotic 
form  in  equation  B19  was  observed  to  approach  the  calculated  val¬ 
ues  for  large  n  and  for  values  of  <|/0  used  for  the  graphs. 

Now  turn  to  the  Legendre  function.  As  v  +  ®  for  fixed  y  and 
0  <  5  <  7i  we  have5B 


P^(cosU)  ) 


r  (  v+m+1)  /n 

r(v  *  if 


sinU) 


) 


-1/2 

cos 


J  +  Hlj  +  o,v-l) 

(B21) 


Thus  as  n  -*■  00  and  remembering  that  n  is  odd  we  have 


p>> 


Tfrff)1'2^  [(n  ♦ 11  *]  ♦  0(n'1) 


n+1 


=  (-1) 


(2 \1/2  r(n  + 

w  rjn  +  fj 


I \r  +  0  (n**1) 


(B22) 


The  gamma  functions  have  an  asymptotic  form  known  as  the  Stirling 
approximation's  which  has  as  v  ■+  00 


5b. 

Ref. 

tt 

eqn. 

00 

• 

»-* 

0 

• 

• 

6B . 

Ref . 

6, 

eqn. 

6.1.37. 
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r(v)  =  e"V  v  2  (2tt)  1/2  [1  +  0  (v  l)] 


(B23 ) 


Applying  this  to  the  gamma  functions  in  equation  B22  gives  as 


n 


T(n  +  2) 


n+y 

- - ^n  +  21 - [l  +  oln'1)] 

e^L  +  *rl 


r("!)  .-+(..§) 


-i  (n+|)ln{n+2)  - (n+1) In (n+|)  . 

=  e  e'  v  e  \  2/  [1  +  0(n-1)] 


=  exp  { (n+1) In 


=  exp  (n+1) 


ln/^±|\  +  \  In  (n+2 )  -  ~  [1  +  0(n_1)] 

H/  I 

+  o (n~2)  +  |  In (n+2 )  -  j  [1  +  0(n-1)] 

ln+l  J 


=  exp | In t (n+2) 1/2 ]  +  0(n_1)|[l  +  0(n-1)] 


=  (n+2)1/2  e0(n  }  [1  +  OCn'1)] 


=  n1/2[l  +  0 (n-1) ] 


Thus  for  the  Legendre  function  as  n  ■>  ®  we  have 


n+1  , /0 

P„(0)  =  (-D  2  n1/2[l  +  0  (n-1)  ] 


[p>>]2  =  |  n[l  +  0  (n-1)  ] 


(B24 ) 


(B25) 
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Next  are  the  Bessel  functions.  As  n  <»  we  have 


kaj  (ka)  ..  .  n+1  .  ,  ..  .  n  ,  -1 

n  _  (ka)  ri  ,  rt/„~l\i  (n+1)  (ka)  ri  ,  nt~.~ l\i 

[kaj  (ka)  ]  1  “  (2n+l)  i  I  11  0(n  )]  “(2"n+l)  !  !~[1  +  0  (n  )] 


HTr[1  +  0{n_1)] 


=  ~[1  +  0  (n-1)  ] 


(B26) 


kah^2) (ka) 


r.  .  C 

Lkahn 


+  °(n~1)]  (“n>  l(2n~nT~r  [1  + 

(ka) n  |  (ka)n+1 


j^[l  +  0(n ~1)] 


(B27) 


Now  we  can  consider  the  terms  in  the  sums  for  the  normalized 
admittances.  As  n  +  “>  we  have 

y<„.  -  i*  I  "ii  +  I  ti  +  om"1)) 


No]' 


[  TTNlpo  ] 


larccos 


(i-2^)  -  j]  +  °(n  2) 


-  1  +  0(n'1)1[1  -  *o]1'/2[iiir]1/2 

•  jcos^j  +  ^arccos^l  -  2^2^  -  +  0(n  1) 
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(B28) 


This  last  result  is  also  the  asymptotic  form  for  Yextn  f°r  large 
n.  Thus  the  individual  terms  in  the  sums  fall  off  like  n~3/2  for 
large  n. 

Equation  B8  gives  the  truncation  error  in  stopping  the  ad¬ 
mittance  sums  at  n  =  ni.  Define  an  upper  bound  for  these  errors 
by 


*1  >  max(|Aint|,  |Aext|> 


(B29) 


Using  the  asymptotic  form  developed  in  equation  B28  we  can  give  an 
approximate  value  for  A^  by  setting  the  cosine  to  one  and  summing 
the  magnitude  of  the  dominant  term  giving 


T  2l1/2r  2  11/2  °°V^ 

*1  *  M1  -  *2]  14:]  J>2 


(B30) 


Replacing  the  sum  by  an  integral  we  have 

3 


r  2l1/2r  2  I1/2  1  C 

ai  ■  m1  -  *2]  k{ 


CO  — 

2  , 

v  dv 


r  21 1/2  r  2 i1//2  ~i 

=  4kal1  -  *oJ  [4;]  " 


r  2i1/2r  2  i1//2 

=  4kati  •  ♦»]  fe] 


which  for  small  iJj0  is 


A^  =  4ka 


[**2n2] 


1/2 


(B31) 


(B32) 


Actually  this  is  a  rather  large  overestimate  of  the  truncation  er¬ 
ror  since  the  cosine  function  alternates  in  sign  and  it  was  re¬ 
placed  by  one. 
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Appendix  C:  Numerical  Techniques  for  Computer  Calculations 
Joe  P.  Martinez,  Dikewood  "  — 

The  numerical  results  plotted  in  figures  4  and  6  through  14 
were  calculated  using  the  Control  Data  Corporation  6600  Computer 
in  the  Air  Force  Weapons  Laboratory  at  Kirtland  AFB. 

The  calculations  for  the  short-circuit  current  transfer  func¬ 
tion  and  the  admittances  included  the  use  of  the  spherical  Bessel 
functions.  These  functions  were  computed  by  using  both  forward 
and  backward  recurrence  techniques.  During  the  initial  attempts 
at  producing  the  Bessel  functions  of  the  first  kind,  the  jn(ka), 
machine  round-off  error  was  encountered  with  the  use  of  the  for¬ 
ward  recurrence  relationship  as  well  as  by  series  representation. 
Backward  recurrence  was  found  to  be  one  solution  in  overcoming 
this  problem.  In  using  this  technique,  where  N  is  the  largest 
order  desired,  the  (N  +  6)th  function  is  set  equal  to  zero,  the 
(N  +  5)th  function  is  set  equal  to  10~50/  and  the  recurrence  rela¬ 
tionship, 


Wka)  “  TT1  3n(ka>  -  Wka>  (cl 

is  then  used  to  determine  the  functions  to  n  =  1.  A  ratio  of  the 
exact  value  to  the  calculated  value  at  n  =  1  is  taken  and  all  the 
calculated  values  to  n  =  N  are  multiplied  by  this  ratio.  The 
Bessel  functions  of  the  second  kind,  yn(ka),  are  calculated  by 
finding  the  first  two  orders  and  then  applying  the  recurrence  re¬ 
lationship 


yn+l(ka)  =  yn(ka>  '  yn-llka)  (C2) 

By  comparing  the  values  produced  by  these  methods  against 
tables  found  in  handbooks,  it  was  determined  that  accuracy  to  eight 
significant  figures  was  obtained  over  the  range  of  n  and  ka  used  in 
the  computations  for  this  note. 

Legendre  functions  of  order  1  are  used  in  the  calculations  of 
equation  72.  These  were  calculated  by  computing  the  functions  of 
degree  1  and  2  and  then  applying  the  recurrence  formula 

PLlM  =  (2  +  s)  X  pn(x)  -  (2  +  k)Pn-lix>  (C3) 

where  x  =  cos(8^).  Results  obtained  by  this  method  were  compared 
with  handbook  tables  and  with  calculations  done  by  series  repre¬ 
sentations  and  backward  recurrence.  These  values  were  found  to  be 
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accurate  to  eight  significant  figures  for  the  range  of  n  used  in 
the  computations. 

While  performing  the  calculations  for  the  short-circuit  cur¬ 
rent  transfer  function  (equation  72)  it  was  determined  that  the 
series  was  a  rapidly  converging  one  and  that  it  met  the  ratio  test 
for  convergence.  This  meant  that  an  upper  bound  on  the  absolute 
error  could  be  expressed  as^C 
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(C4 ) 


where  aft  is  the  last  term  included  in  the  summation.  As  explained 
in  Section  IV,  only  the  odd-numbered  n  are  included  in  the  summa¬ 
tion.  The  relative  error  will  then  be 
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In  the  calculations,  E  was  set  at  .001.  The  number  of  terms  nec¬ 
essary  to  satisfy  this  criterion  varied  with  ka,  the  larger  values 
requiring  more  terms.  Generally  the  error  criterion  was  satisfied 
in  fewer  than  100  orders  (50  terms) . 

The  calculation  of  the  yint  and  yext  sums  (equations  105  and 
110)  were  handled  in  a  different  manner  than  that  of  the  T  sum  due 
to  their  slow  convergence.  The  behavior  of  the  expressions  involv¬ 
ing  the  spherical  Bessel  functions  and  the  hypergeometric  function 
were  studied  separately  to  determine  how  the  numerical  calculations 
could  be  made. 

It  was  shown  in  appendix  B  that,  as  n  +  ® 


kaVka)  ka 

[kajn (ka) ] '  n  +  1 


and 


kah  ^  (ka)  . 

[kah<2>(ka)]'  * 


(C6) 
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1C.  W.  Kaplan,  Advanced  Calculus,  Addison-Wesley ,  1952,  pp.  328, 
329. 
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Numerical  calculations  were  made  in  order  to  determine  at  what 
point  the  asymptotic  forms  of  these  expressions  could  be  used  with 
minimum  loss  of  accuracy.  It  was  desirable  to  switch  to  the  as¬ 
ymptotic  form  at  some  point  to  avoid  round-off  numerical  errors 
due  to  "he  extremely  small  magnitude  of  the  jn(ka)  at  large  n.  It 
was  found  that  at  n  =  10|ka|  +  30  the  asymptotic  forms  may  be  used, 
with  a  relative  error  of  .03%  at  ka  =  1.0  introduced  in  the  terms. 
Of  course,  as  n  gets  larger  this  error  will  diminish  with  each 
term. 


In  performing  the  numerical  calculations  for  the  hypergeo¬ 
metric  function  of  equations  105  and  110  it  was  found  that  numer¬ 
ical  round-off  error  was  introduced  after  the  sum  of  equation  103 
was  carried  to  a  large  number  of  terms,  tne  number  depending  on 
\JjQ.  Going  to  double  precision  on  the  computer  greatly  improved 
the  situation,  as  now  the  number  of  usable  values  was  doubled. 

But,  the  round-off  errors  were  still  introduced  before  the  yint 
and  yext  sums  could  be  completed.  As  was  mentioned  in  appendix  B, 
the  asymptotic  form  (equation  B19)  approaches  Fn  for  large  n.  So 
it  was  decided  to  switch  to  this  form  just  before  the  numerical 
errors  started  to  be  significant.  Table  Cl  is  a  summary  of  the 
pertinent  values  which  occur  immediately  before  the  switch-over 
points. 

Relative 
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\  Asymptotic 
}  Formula 

Absolute 

Difference 

Difference 

(Percent) 

.001 

24999 

. 055569 

.055816 

.000246 

.443453 

.01 

3199 

.093564 

.093394 

.000170 

.182042 

.  1 

319 

.095720 

.096555 

.000835 

.872126 

Table  Cl.  Values  of  the  hypergeometric  function 
one  term  before  the  switch-over  point 


Now  that  methods  have  been  established  to  calculate  the  terms 
involving  the  Bessel  functions  and  the  hypergeometric  function, 
the  calculation  of  the  yint  and  Yext  sums  may  proceed.  At  some  ng 
a  switch  to  the  asymptotic  forms  of  equations  C 6  and  C7  takes 
place,  so  that  equation  105  becomes 
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Similarly,  equation  110  becomes 
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Since  the  switch-over  occurs  when  n  >_  10|ka|  +30,  and  the  largest 
value  of  ka  used  in  these  calculations  was  10.0,  all  changeovers 
to  the  asymptotic  forms  will  take  place  by  n  >  131.  It  was  found 
to  be  convenient  to  calculate  the  sums,  with  Iffka  factored  out, 
from  n  =  201  to  00  (n  odd)  first,  and  then  include  these  values  in 
the  calculations  later.  Since  ka  can  be  factored  out  from  the 
last  part  of  the  summation,  these  values  may  be  used  with  any  ka. 
The  next  problem  was  that  of  determining  how  far  n  should  be  car¬ 
ried.  This  was  done  numerically  by  running  the  sums  out  progres¬ 
sively  further  to  larger  n  and  noting  the  relative  change.  The 
sums  for  all  ipQ  were  carried  to  n  =  100001  and  it  was  noted  that 
they  were  still  changing  about  1  part  in  1000,  or  .1%.  It  was 
further  noted  that  there  were  peaks  and  minimums  occurring  as  the 
sum  progressed,  forming  a  diminishing  envelope  converging  at  some 
number.  The  last  few  maxima  and  minima  were  then  taken  and  the 
value  to  which  they  were  converging  was  extrapolated  to  5  signifi¬ 
cant  digits.  Table  C2  gives  these  values.  The  relative  error  in 
the  sum  is  then  in  the  order  of  .05%.  This  error  is,  of  course, 
multiplied  by  irka  when  the  sums  are  used  in  the  admittance 
calculations . 
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ext 
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.001  1.0987  x  10° 

.01  -8.6355  x  10"2 

.1  -2.0147  x  10-3 


1.1013  x  10° 
-8.6605  x  10-2 
-2.0234  x  10-3 


Table  C2 .  Values  to  which  admittance  sums  converge 
from  n  =  201,  with  iffka  factored  out 


Most  of  the  other  numerical  calculations  performed  for  the 
graphs  in  this  note  are  straight  forward  and  no  explanation  is 
required . 
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